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Abstract 

We investigate whether known properties of the rj' meson are consistent with its 
being the Goldstone boson of the spontaneously broken anomaly-free R symmetry 
required in the light gluino scenario. We fit the masses and 27 decays of the r/ and 7/ 
mesons, and also their production in radiative J/ ip decays. We find that the 77 — 7/ 
system is well-described in the light gluino scenario, if ni\ ~ (84 — 144) MeV and 
(AA) ~ -(0.15-0.36) GeV 3 . These values are in the range expected when the gluino 
gets its mass entirely from radiative corrections. 
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It has been known for a long time that there exists a window for a long-lived 
light gluino in SUSY phenomenology |l], 0. Dimension-3 SUSY breaking operators 
(including gaugino masses) may naturally be highly suppressed in the MSSM Q, 
with a number of attractive consequences such as eliminating the CP problem of the 
MSSM J?], |[] and providing the observed dark matter ||. Gauginos acquire masses 
through radiative corrections from electroweak and top-stop loops || [1. Estimates 
of the constrained soft SUSY breaking parameter space lead to the gluino mass range 
nix — (0.1 — 1) GeV |§. Such a light gluino leads to a slower running of the strong 
coupling constant and might help to resolve the strong coupling constant discrepancy 
problem (see || and refs. therein). 

The Lagrangian of supersymmetric QCD after decoupling heavy squark modes 
consists of the usual QCD plus light gluinos. Hence, among the usual symmetries 
of QCD there exists an additional chiral symmetry associated with the chiral U(l) 
transformation of the gluino field. A particular linear combination of the currents as- 
sociated with the quark and gluino chiral rotations is free of the gluon anomaly. We 
will call this anomaly- free symmetry R invariance, and the corresponding current the 
R current. The R invariance is spontaneously broken by quark and gluino conden- 
sates. It was argued in ref. JK| that supersymmetry must be explicitly broken by a 
gluino mass of order 1 GeV or more, in order to avoid an unobserved Goldstone boson. 
However it was suggested in ref. JlJ that the required pseudogoldstone boson might 
be identified with the rj'; a naive estimate for the required gluino condensate was 
given in ref. H]. In this scenario the h state which gets its mass from the anomaly 
is expected to have mass approximately 1.5 GeV and is naturally identified as the 
otherwise-mysterious "extra" singlet pseudoscalar at 1410 MeV [|TTJ. Approximately 
degenerate with it is a spin- 1/2 "glueballino" whose detection is discussed in ref. fill . 

Here we investigate in greater detail the possibility that the 7/ is the required 
pseudogoldstone boson of spontaneously broken R-invariance. In addition to giving 
a more accurate determination of the gluino condensate, we examine whether a non- 
vanishing gluino condensate is consistent with the 77 decay rates of 77 and 77', and 
with the relative production rates of 77 and 77' in J ftp radiative decays. As in the 
theory without gluinos, the 77, rj' system must be analyzed together, since mixing 
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between flavor octet and singlet components is essential for understanding their 77 
widths. A basic tool in our consideration is the QCD sum rule method [13[] (reviewed 



in refs. []14| , |15|| , [|T( 

The paper is organized as follows. We begin by fixing the conventions for the 
anomaly free axial currents and mixing angles. After this, we derive expressions 
for the singlet and octet pseudoscalar meson decay constants using the sum rule 
technique. Next we consider the 27 decays of rj and 7/ mesons and derive sum rules 
for J/ip decays into n{ r q') r y. Finally, imposing all available experimental restrictions 
on these processes, we numerically solve our system of equations, leading to estimates 
for the light gluino mass and condensate. 

Let us start with the definition of the anomaly free axial currents 

jjS = 7^ (m7 m 75« + hulsd - 257^755) (1) 

J$t> = 7^(w7m75W + dl^d + s7 M 7 5 s - ^7^5 A), (2) 

where the Majorana spinor A denotes the gluino field. The derivatives of these currents 
play the role of interpolating operators for pure octet and singlet pseudoscalar meson 
states. The physical n and 7/ mesons are defined through the mixing angle 9 for 
which we choose the standard parametrization \rj) = \rj s )cos8 — \i]°)sin8, \n') = 



\rp)cos6 + \r] )sin6. Following ref. ||17| , define the decay constants Fq and Fg as 



(0|<V,2V> = ^0 cosd m\, (0|<v£°5k> = ~*b smO 
r(8)u'\ _ v „A„n ™2 /n , A T (8), 



m 



'/■ 



<o|<v;f 5 v> = ^8 sine mj, (o|<vW = F * cosd K 

Consider now the correlator of two axial currents 

tA,,, = v? J e^(0\TJ^(x)J^(0)\0)d 4 x = -q u Ii{q 2 = -Q 2 ). (3) 
In accordance with the QCD sum rule approach [[13] the Borel transforms of the 



phenomenological and theoretical parts of the correlator are equated, making the 
standard decomposition 

ImU phen (s) = ImU poles (s) + ImU pert {s)e(s - s ), 
ImIl theor '(s) = ImW ert (s) +ImU cond (s). 
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The superscripts "poles" , "pert" , and "cond" label the resonance, perturbative and 
power (condensate) contributions to the correlator. The result is 

/* OO /*5Q P OO 

/ e^ImU poles (s)ds = ImU pert (s)ds + ImU cond (s)ds. (4) 

Jo Jo Jo 

The conventional parametrization for the pole and the condensate contributions are 
the following: 



oo oo 



ImU poles (s) =<ixY J C i 5{s-ml), ImU cond (s) = n ]T kj5^{s)/j\, (5) 

i=l j=0 

with Cj's being pole residues, k/s being gauge invariant condensates and denoting 
the j'th derivatives of the Dirac's delta function. Substituting these expressions into 
eq. (£|) and expanding in inverse powers of the Borel parameter (for M 2 » sq) we 
get the sum rules in each order of 1/M 2 . For practical calculation it is enough to 
keep only the first two equations of this tower. These equations look like 

7r(ci + c 2 ) = f ° ImW er \s)ds + 7rk , (6) 
Jo 

c x m\ + c 2 m\) = / sImU pert (s)ds - irh. (7) 
Jo 



7TI 



It is not difficult to see that these relations are nothing but the QCD finite energy 
sum rule (FESR) [JTSJ , modified by the condensate contributions |Tj| . 

Having set up the general framework, let us turn to the calculations for our ap- 
plication. In the case at hand, the one loop calculation for the singlet correlator in 
eq. (4) leads to the following resultsQ 



v ZTT SIX ' 
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ko = - 2msi ' S) - h = -(^ + ml)(-Gl). 

3 6 v 6 7T 

The constants Ci, c% , mi and m,2 entering the phenomenological part of the sum rules 
are: C\ = Fq m 2 sin 2 9, m\ = m 2 ; c% = Fq m 2 , cos 2 9, m\ = m 2 ,. Substituting all 
these relations into eqs. (|6|) and (0), determining then so from eq. (0), and plugging 

*We hereafter neglect the u and d quark masses. 
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back its value into eq. 
constant: 



), we arrive at the following expression for the singlet axial 



' m 2 s m\ \ (1 + 0.1 tan 2 6)(l + tan 2 
2^ + 3^ J (1 + 0.33 tan 2 6) 2 



x 



1 + 

V N 



2m s (ss)/3 + m x (XX)/6 (1 + 0.33 tan 2 6) 



in 



/m| , m|s (1 



0.1 tan 2 ( 



( 



2m s (ss) 
3m 2 , 



m\{\\) \ (l + tan 2 #) 
6m 2 , / (1 + 0.33 tan 2 0)' 



where 







;^ + ^)<^y (1 + 0.33 tan 2 



4 / < , ™i 



(1 + 0.1 tan 2 #) 2 ' 

Using the results of the above calculation, it is easy to obtain the analogous sum rule 
for the decay constant F$. One just needs to drop the gluino contribution and choose 
the proper normalization for the corresponding axial current. The result is: 



mz 



1 + 9.4 tan 2 6)(l + tan 2 



7T 



'1 + 3.1 tan 2 6) 2 



1 



\ 



1 - 



8vr 2 (ss) (1 + 3.1 tan 2 6) 
3m 2 m s (1 + 9.4 tan 2 6) 



(9) 



-( 



4m s (ss) 



\ 3m 2 



'1 + tan 2 



1 + 3.1 tan 2 ^ 



where 



2n 2 (fG%) (1 + 3.1 tan 2 



3 (1 + 9.4 tan 2 6) 2 ' 

Generally speaking there are five unknowns in these equations. These are the decay 
constants F and F 8 f\ mixing angle 9, gluino mass m\ and the gluino condensate (AA). 
Our strategy hereafter will be to impose the restrictions on these unknowns coming 
from the 2j decays of the n and rf mesons and from the radiative decays of the J/ip. 

Let us turn first to the 27 decays of the n and n' mesons. These decays are 
governed by the electromagnetic axial anomaly. Since gluinos do not interact with 
photons in leading order, the gluino part of the singlet current just drops out from the 
calculation of the corresponding triangle graphs. Therefore the formal expressions for 



t Though Fg is known from chiral perturbation theory pi[ , we include it in the list of unknowns 
to provide a consistency check of our calculations. 
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the decay widths are the same as in conventional QCD. However the interpolating 
current for the rf meson is now free of the gluon anomaly and one need not worry 
about modification of PCAC as is necessary in QCD fl22| . (Because this current is 
already modified by subtracting the proper amount of gluinos). Taking into account 
mixings, one obtains 



r(?7->-27) _ /mgw F^cogg 2F^sin0 vi 
r(7r^2 7 )~ Vm^iW6 F Vs > 1 J 

T(r]' -> 27) _ /T^w ggroogg F^sinO ^i 

r( 7 r-2 7 ) ^mjA Fq ^3 + F 8 VQ ] ' 1 J 

The available experimental data for these decay rates are ||23|| : T(rj' — > 27) = (4.26 ± 
0.62) keV, T(r] -> 2 7 ) = (0.51 ± 0.09) keV, T(7r -> 2 7 ) = (7.74 ± 0.58) eV. 

Eqs. ( |T0D and ( |TT|) are derived using PCAC and soft meson technique. The latter 
is a dubious approximation for the rf mesonf]. To obtain some estimate of the effect 
of the extrapolation from m v r = to m v > = 958 MeV, we apply the interpolation 
technique of ref. |[24| . The result is that the Fo's appearing in eq. (|Tl|) are multiplied 
by the factor (1 — (m^//167r 2 F 2 )). Results of both estimates will be reported below. 

Now let us turn to the radiative decays of vector quarkonium. In particular we will 
concentrate on J/ip — ► 77 (7/) 7. The basic properties of these processes were worked 



out using the sum rule technique in refs. pq| , [ 26[ . To understand the mechanism 
which dominates the decay, it is relevant to look at the quantum numbers of the 
particles involved: J/ip{l ) — > PS(0~ + ) 7(l p ~). In general the spatial parity of the 
emitted photon may be negative (for "electric" transition) or positive (for "magnetic" 
transition - see for example [f27j). From the quantum numbers of the quarkonium and 
the pseudoscalar meson, the emitted photon in the present case must have positive 
spatial parity and consequently should be the "magnetic" type. Keeping this fact in 
mind, the following three step picture for these decays emerges naturally: the heavy 
cc quark system being in a vector state emits the "magnetic" photon and turns into 
a pseudoscalar state of two c quarks (a virtual r] c meson). Then this virtual r] c emits 
two gluons in a pseudoscalar state which hadronize into the final pseudoscalar meson 
(77 or rj'). Using the operator ic^^c for the interpolating field of the r] c meson, one 
* We are grateful to H. Georgi for bringing this point to our attention. 
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obtains the following ratio: 



r(J/ip -> 7/7) (0\icry 5 c\rf) 



2 /vr^T,, — ml, S 3 



r(J/V^?77) (0\ic<y 5 c\ri) 

It is convenient to introduce the notation (0|ic75d|?/) = a 2 a v //iT 2 and (0|zc75c|r]) = 
a 2 a v /7i 2 . Eq. (|12|) and the experimental data on these decays [^| give 

(2.48 ±0.16)|cg. (13) 



\a v >\ 



The next step is to relate a v and a'^ to properties of the r\ and rf using sum rules. 
Therefore consider the two-point correlator 



P( g ^) = i y e^{0\Tiey 5 c(x)d^(0)\0)d*x. (14) 

The calculational procedure is the same one used above for the determination of Fq 
and F 8 . The phenomenological part is saturated by the physical 77 and 7/ meson 
states. The leading perturbative part is determined by the three-loop diagram shown 



in fig. 1(a). This diagram has been calculated recently in ref. [28], using an expansion 



in inverse powers of heavy quark mass [ 29] . In our case this is the leading approxima- 
tion in l/m c . The leading contribution to the nonperturbative part arises from the 
diagrams of fig. l(b,c,d). We have used the heavy quark mass expansion and dimen- 
sional regularization, with the 't Hooft-Veltman prescription |3(| for the 75 matrix 
in D = 4 — 2e dimensional space-time. The calculational technique is based on the 



integration by parts method of ref. [ 3 1 



The results corresponding to the three diagrams of figs, lb, lc and Id are 
Im P b (s + is) = 0; Im P c (s + ie) = Im Pd(s + is) = 

1 :( ^ 2 ( !^_l!!^v (2G(1)1) _ G(2)1))) 



2\/3 ^ vr ' 3m c 24 m c 
where G(l, 1) = 1/e + InAir - 7 + 2 + o(e) and G(2, 1) = -1/e - InAir + 7 + o(e) are 

the G-functions introduced in ref. and 7 ~ 0.5772 is the Euler constant. Using 

these results for the theoretical part, the sum rules are 

F (a v/ m 2 f cos6 - a v m 2 sin6) = — jfi- 3 f f{s)ds - 

m c \/3 87r z Jo 

1 m s (ss) 9 m A (AA) A 

7! ( "^r + 8^^) So ' (15) 




/HK yHK ATx 

XQQO/' xQQO/' xfiQO/' 



Figure 1: The double plain lines denote the c quark propagators and the single plain 
lines denote the light quark or gluino propagators. 



Fo(a r) /mtcos9 — a^mtsinO) 



2(m 2 -3ml) 3 r a ° 



m, 



.VS 8tt 2 Jo 



sf(s)ds 



J_ _ m s (ss) 9 m x {W) \sl 
m c 8 m c ' 2 : 



(16) 



where the function f(s) = As — slnA^ + ^ (^ — Yos^^f) appears as the imaginary 
part of diagram (la) p8| . 

In analogy with the singlet case one can write down the same sum rules for the 
correlator with the nonsinglet current. Taking the proper normalization for this last 
we get 



Fg(a v m 2 cosO + a v iin 2 ,sint 



- im 2 s 3 [4 8) 

1 f(s)ds - 



.V6 8n 2 Jo 



1 2m s (ss) ( 8) 



(17) 



Am 2 3 r'o 



(8) 



F 8 (a v mcos8 + a v im,sin8) = p— — / sf(s)ds — 

1 1 m c v6 o7r z Jo 

1 2m s (ss)s { m 

rn c 2 1 ' 

We now have enough equations to determine the unknown constants. In fact, we 



have nine equations (|-[TT|), fll3|) , (p~5| - p~8D and nine unknowns (6, F , F 8 , m\, (XX) 



cy, So and s^) which depend on m s , (ss), (^-G 2 ). For the s quark mass, recent 
lattice calculation gives m s (2 GeV) = (141 ± 17) MeV 1)331 ; QCD sum rules give 
m s (l GeV) = (189 ± 32) MeV @ and m s (l GeV) = (171 ± 15) MeV g|. In our 
calculations we take the average of the above three results m s (l GeV) = (167 ± 13) 
MeV. We take the s quark condensate from ref. |20[ (see also the recent paper f36| ) 
(ss) = (0.7±0.2)(uu) = -(0.011 ±0.003) GeV 3 . There are a number of estimates for 



the gluon condensate in the literature (see refs. [p^,fll3|). We take the world average 
value of these calculations (fG 2 ^) = (2.5 ± 0.9) 10~ 2 GeV 4 . 

Our procedure is this. First we use eq. (UJ to find F (or F (l - (ml, /16tt 2 F 2 )) 
using the interpolation based on the Brodsky-Lepage analysis) and substitute it into 
eq. (p"l~D. After this we have two equations, (|9|) and (|Tl|), for F$ and 9. Solving this 
system numerically we find F 8 ~ (130±20(stat.)) MeV and 9 ~ -(23.5° ±2.7° (stat.)). 
The value obtained here for F 8 agrees within error bars with the chiral perturbation 
theory estimate F 8 ~ 1.25F n ~ 116 MeV |[21|| . Our value for 9 is also in agreement 
with previous calculations (see for example ref. [0). This gives confidence in these 
values for m s and (ss), and in the present method combining PCAC and sum rule 
techniques for calculating the decay constants. We find F ~ (93 ± 17(stat.)) MeV 
with the conventional soft meson extrapolation for the //, and Fq ~ (136 ± 13(stat.)) 
MeV with the alternative procedure described above. 

Next we make use of the data and SR's for J/ip — > 7/7 and J/ip 7/7 , namely eqs. 
(|l3| , p!7| , [T8|) to determine a Vl a v > and numerically using Mathematica. The proce- 



dure is trivial but cumbersome and yields: \a v /\ ~ (202.1 ± 45.3)10 3 GeV 2 , \a v \ ~ 
(81.6 ± 18.1)10- 3 GeV 2 and s 8) - (1-7 ± 0.2)GeV 2 . Note that we took from experi- 
ment the ratio laj/lflyl, but have not obtained their individual values. Thus if we 
had a reliable independent way of estimating the matrix element (0|zc7sc|?7 c ) = F c 
appearing in the amplitude for J/ip — > 7/7 and J/ip —>■ rj'^, we would have a predic- 
tion for the overall magnitude of these rates. Unfortunately, F c cannot be determined 
independently with sufficient accuracy to give us a useful testQ. 

The results above allow us to turn to the calculation of the basic quantities, the 



§M. Shifman, private communication. 
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gluino mass and condensate. Substituting the quantities obtained earlier into eqs. 
(H [15], [161), we have a system of three equations with the last three unknowns: 
m\, (XX) and sq. The numerical solution of this set provides us with the fol- 
lowing results using Fq extracted with the conventional soft meson aproximation: 
s ~ (1.8 ± 0.2) GeV 2 , m A ~ (120 ± 24) MeV, and (AA) ~ -(0.22 ± 0.07) GeV 3 . 
Using F obtained from the alternative extrapolation in m v i gives: So — (1.9 ± 
0.2) GeV 2 , m x ~ (102 ± 18) MeV, and (AA) ~ -(0.31 ± 0.05) GeV 3 . From this 
range of results we arrive at the final estimate 

m A ~ (84 - 144)MeV, (AA) ~ -(0.15 - 0.36) GeV 3 . 

Before concluding, we note that in conventional QCD the rj' is identified with the 
particle which acquires mass due to the axial anomaly. The Witten-Veneziano (WV) 
formula |41| relates the vacuum topological susceptibility of pure Yang-Mills theory 



to the decay constant and mass of the rf meson in the large N c limit. (For a recent 
discussion see ref. Ift^l). In the limit when at least one quark is massless this relation 
looks like 

/ 2 \ 2 1 



J (0\TG a G a (x)G a G a (0)\0)d 4 x\] v ^ oo = -F^U^. (21) 



V167T 2 , 

In the scenario we have considered here, the 7/ is a pseudogoldstone boson and a 
different state acquires its mass from the anomaly. (In supersymmetric YM theory 
this state is a pure gluino-gluino bound state j|t|.) It is possible to derive an analogous 
formula to eq. (21) in the theory with a light gluino. The only difference is that the 
left hand side should be calculated in supersymmetric YM theory. The right hand 
side is unchanged, and still depends on the decay constant and mass of the physical 
7/ meson which now contains a light gluino. Unfortunately, the accuracy of the 
determination of the vacuum topological susceptibility (for a review see ref. ||42|| ) is 
not sufficient at present to decide which of the two scenarios is realized. A detailed 
discussion of these and related topics will be presented elsewhere. 

To summarize, in this paper we have undertaken an attempt to identify the 7/ 
meson with the Goldstone boson of the spontaneously broken, anomaly-free R sym- 
metry arising in a theory with a light gluino. We have calculated the corresponding 
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decay constant Fq using the QCD sum rule method. As a consistency check, the 
analogous calculation for the octet decay constant F s is performed. The result is in 
good agreement with what was already known before. The experimental values for 
the decay rates of r]'(j]) — > 27 and J/ip — > 7/(77)7 impose certain restrictions on the 
values of gluino mass and condensate. Using the whole system of equations coming 
from sum rules and from quantities determined by experiment, we determined the 
gluino mass and condensate. The value of the gluino condensate indicates that the 
scale of the breakdown of /^-symmetry is about 2.5 times larger than that of chiral 
symmetry. In view of the factor 9/4 larger value of the adjoint Casimir compared 
to the fundamental one, this seems to be of the right magnitude. Note also that the 
gluino mass value we have found, m\ ~ 80 — 160 MeV, is in the range predicted in 
ref. 0. If the properties of the 77, 77' mesons had been incompatible with the light 
gluino scenario, our procedure would have given nonsensical rather than reasonable 
values of gluino mass and condensate. Thus our result shows a surprising consistency 
between the properties of the ground-state pseudoscalar mesons and the light gluino 
scenario, and favors a gluino mass of order 100 MeV. 

The authors are grateful to K.G. Chetyrkin for discussion of the results of ref. 



8| and M. Losada for help in using the "feynMF" package p9 
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